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We describe a simple dynamical model characterized by the presence of two noncommuting Hamil- 
tonian constraints. This feature mimics the constraint structure of general relativity, where there is 
one Hamiltonian constraint associated with each space point. We solve the classical and quantum 
dynamics of the model, which turns out to be governed by an SL(2, R) gauge symmetry, local in 
time. In classical theory, we solve the equations of motion, find a SO(2, 2) algebra of Dirac ob- 
servables, find the gauge transformations for the Lagrangian and canonical variables and for the 
Lagrange multipliers. In quantum theory, we find the physical states, the quantum observables, and 
the physical inner product, which is determined by the reality conditions. In addition, we construct 
the classical and quantum evolving constants of the system. The model illustrates how to describe 
physical gauge-invariant relative evolution when coordinate time evolution is a gauge. 
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I. INTRODUCTION 

General relativity (GR) has a characteristic gauge in- 
variance, which implies that its canonical Hamiltonian 
vanishes weakly. As a consequence, its dynamics is not 
governed by a genuine Hamiltonian, but rather by a 
"Hamiltonian constraint" . This peculiar feature of the 
theory has a crucial physical significance, connected to 
the relational nature of the general-relativistic spatiotem- 
poral notions and raises a number of important 

conceptual as well as technical problems, particularly in 
relation to the quantization of the theory m. In the past, 
much clarity has been shed on these problems by study- 
ing finite dimensional models mimicking the constraint 
structure of the theory, and in particular, having weakly 
vanishing Hamiltonian ||. 

There is an aspect of the constraint structure of GR, 
however, which, as far as we are aware, has not been an- 
alyzed with the use of constrained models. In GR, there 
isn't just a single Hamiltonian constraint, but rather a 
family of Hamiltonian constraints, one, so to say, for each 
coordinate-space point. Furthermore, the Hamiltonian 
constraints do not commute with each other (have non- 
vanishing Poisson brackets with each other). Indeed, the 
constraint algebra of GR has the well known structure 

{H, H} ~ D, {H, D} ~ H, {D, D} ~ D, (1) 

where H represents the Hamiltonian constraints and D 
the diffcomorphism constraints. In this paper we present 
a model that mimics this aspect of GR. 

The model we present has three constraints, which we 
call Hi, H2 and D. Their algebra has the structure 



{Hx,H 2 }r.D, {H Z ,D}~H Z 



(2) 



which mimics ([!]). (Models with several commuting 
Hamiltonian constraints were considered in §.) The con- 
straints H\ and H2 are quadratic in the momenta, while 
D is linear, as their correspondents in GR. 

The model has an interesting structure which exem- 
plifies in a non trivial manner various aspects of the 
quantization and interpretation of the fully constrained 
systems. We analyze in detail its classical and quan- 
tum dynamics, which can both be solved completely. We 
display the general solution of the equations of motion 
and the finite gauge transformation of variables and La- 
grange multipliers. The constraint algebra turns out to 
be SL(2, R) and the model is invariant under an SL(2, R) 
gauge invariance, local in time. We find a complete 
50(2, 2) algebra of gauge invariant observables, as well as 
a (smaller) complete set of independent observables. The 
phase space turns out to have the topology of four cones 
connected at their vertices. We then study the quantum 
dynamics, solve the Dirac constraints, exhibit the phys- 
ical states explicitly, and construct a complete family of 
gauge invariant operators. The reality properties of the 
gauge invariant operators fix uniquely the physical scalar 
product. In addition, we define the classical and quan- 
tum evolving constants |6) of the system, and we discuss 
the observability of evolution for the systems (like GR) 
in which time is a gauge and the theory has no preferred 
physical time. 



II. CLASSICAL DYNAMICS 



Definition of the model. The model we consider is de- 
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fined by the action 
S[u, v, N, M, A] = 

dt [ N (Vu 2 + v 2 ) + M (Vv 2 + u 2 ) ' 



where 



Vu 



— (u-Xu), 



Vv 



— + Xv) ; 



(3) 



(4) 



the two Lagrangian dynamical variables u — (it 1 , u 2 ) and 
v = (w 1 ,^ 2 ) are two-dimensional real vectors; N, M and 
A are Lagrange multipliers. The squares are taken in R 2 : 
u 2 = u ■ u = (u 1 ) 2 + (u 2 ) 2 . 

Hamiltonian analysis. The Hamiltonian analysis is 
simplified by first rewriting the action in the following 
form 



S= dt 



it • Vu + 'v ■ Vv - N^(Vu 2 -v 2 ) 
-m\^Vv 2 - u 2 ) - X(u ■ Vu - v ■ Vv) 



(5) 



From this form, we see that the momenta conjugate to u 
and v are 



p = Vu 



and 



Vv 



(0) 



respectively, and that we have a weakly vanishing Hamil- 
tonian and three primary constraints 



Hi 


1 

~ 2 


if 


- v 2 ) , 


H 2 


1 

~ 2 


(n 2 


-u 2 ) 


D 


= u 


■P- 


- V ■ 7T 



The Hamilton equations of motion are 



u = Np + Xu . 
p= Mu — Xp, 



v = Aln — Xv . 
7T = Nv + Xtt . 



(7) 



(8) 



Using 
straints 



and M) we find the evolution of the con- 



i?i = MD - 2Affi , 
H 2 = -ND + 2XH 2 , 
D = -2MH 2 + 2NH X 



(9) 



These equations show that there are no secondary con- 
straints, and that three constraints (Q) are first class. 
The dynamics of the system is given entirely by the con- 
straints and the Hamiltonian is H — NHi + MH 2 + XD. 
Since we have four real dynamical variables {u and v) 
and three first class constraints, the system has a single 
physical degree of freedom. 

The Poisson algebra of the constraints can be read di- 
rectly from (||); it is given by (cfr. eq. (^)) 



{H x ,H 2 } = D, 
{#! , D} = -2H X , 
{H 2 , D} = 2H 2 . 



(10) 



This algebra is isomorphic to sl(2,R), the Lie algebra of 
the group SL(2,R). 

Analogy with GR. The model has a structure recall- 
ing GR. The analogy is transparent in the Hamiltonian 
framework, given the similar structure of the two con- 
straint algebras. In the Lagrangian framework, com- 
pare the action (|^) with the Einstein-Hilbert action Sqr- 
Written in terms of the Arnowitt-Deser-Misner (ADM) 
variables, Sgr is 



S GR {g,N, A]= fdtf ^fgdx N (Dg 2 + R\g]), 



Dg a b = -j^(gab - 2D (a X b) ) 



(11) 



where g is the three-dimensional metric, N the lapse and 
A the shift, R the three-dimensional Ricci scalar, we have 
indicated the extrinsic curvature by —Dg a }, and written 
Dg 2 = Dg c c Dgj-Dg ab Dg ab . Notice that the two compo- 
nents of u mimic the metric in a space point, the two com- 
ponents of v mimic the metric in a second space point, 
N mimics the lapse in the first point, M the lapse in a 
second point and A the shift. The sum in (^|) mimics the 
integration over x in (|TT|), and the definition of Dv and 
Dv mimics the extrinsic curvature. 

Gauge invariance. Under an infinitesimal gauge trans- 
formation generated by infinitesimal time dependent pa- 
rameters n(t),m(t),l(t), the canonical variables trans- 
forms as [0 

Su = l(t)u + n(t)p, 
dp = m(t)u — l(t)p, 
5fr = l(t)n + n(t)v , 
dv = m{t)'K — l(t)v , 



(12) 



while the Lagrange multipliers transform as 

5N = h(t) - 2n(t)X + 2l(t)N , 
SM = m(t) + 2m(t)X - 2l(t)M , 
SX = l(t) + n(t)M - m(t)N . 



(13) 



We can check the transformation of the action (|5|) under 
this infinitesimal variation of the canonical variables and 
the Lagrange multipliers. We find that 6S — provided 
that the boundary term n(t)(p 2 + v 2 )+m(t)(n 2 +u 2 )\ t t Z t / i 
vanishes. 

The problem of finding the finite gauge transforma- 
tions can be solved by using the fact that ( |l2| ) is an 
infinitesimal SL(2, R) transformation. More precisely, 
each one of the four pairs (u 1 ,^ 1 ), (u 2 ,p 2 ), (^i; 1 ), 
(7r 2 ,w 2 ) (notice that the order is inverted in the second 
two), transforms in the fundamental representation of 
SX(2, R). It follows that the finite gauge transformation 
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of the canonical variables generated by the first class con- 
straints are given by finite SL(2, R) transformations as 
follows 



u = a(t)u + f3(t)p, 

f = -r(t)u + 6{t)p, 



if' = a(t)if + (3(t)v . 
if = 7(t)7r + S(t)v , 



where the matrix 



G(t) 



_ (a{t) (3(t) 
7 (t) S(t) 



(14) 



(15) 



is in SL(2,R), that is, with the only restriction that 
a(t)S(t) — f3(t)y(t) = 1 . Thus, the system is invariant 
under an SL(2, R) gauge invariance local in time. 

The finite transformation law for the Lagrange multi- 
pliers can be found from the definitions of the momenta. 
We obtain with some algebra 

N' = a 2 N - P 2 M - 2a(3\ + a$ - a/3 , 
M' = -j 2 N + S 2 M + 2jSX + jS - jS , 
A' = -a 7 7V + (38M + {aS + /? 7 )A + a5 - /3j. (16) 

Below we give a clean geometric interpretation of these 
ugly- looking transformations. 

We can now check the invariance of the action. By 
plugging ( |l4| ) and (16) into the action (j|) we get with 
some algebra 



S'= dt 



u-p + il-Tz- (NHi + MH 2 + XD) 



+ 



(/3-y)(u ■ p + v ■ if ) + -(aj)(u 2 + n 2 ) 



v 2 ) 



-\t=tf 



(17) 



The action is invariant provided that the boundary term 
vanishes. 

Solution to the equations of motion. The evolution of 
the system can be viewed geometrically. Let us focus on 
the (u,p) sector -the (if, v) behaves in the same manner. 
The equations of motion (^) for this sector can be written 
in the form 



d 
dl 



X 

M 



N 
-X 



= 0. 



The matrix composed by the Lagrange multipliers is val- 
ued in the Lie Algebra of the SL(2, R) group and can 
be viewed as the Yang-Mills connection for the local (in 
time) gauge group SL(2, R) 



A(t) 



_ ( A(i) N(t) \ 
\M{t) -X(t)J 



(19) 



This is not a vague analogy: using this notation, the ugly 
transformation (flBf) becomes 



(20) 



A' = GAG' 1 - G—G' 1 
dt 



That is, A transforms precisely as a connection. Under a 
time dependent gauge transformation G(t), (u,p) trans- 
form as in ([l4|), A transforms as in (^0|) and the form of 
the equation of motion (|l8|) is preserved. 

Given the geometric analogy, it is easy to integrate 
the equations of motion. The Lagrange multipliers can 
be chosen as arbitrary functions of time, namely we can 
choose an arbitrary time dependent sl(2, R) matrix A(t). 
The solution of the equations of motion (Tlq) is then ob- 
tained from the initial value (uq,Pq) at time t = by 



( u(t) \ _ ( a(t) b(t) \ ( uo 

\p(t) U(t) d(t) \p 



where the matrix 

U(t) 



a(t) b(t) 
c(t) d(t) 



satisfies the parallel transport equation 



-U(t)-A(t)U(t)=0. 



(21) 



(22) 



(23) 



The solution is given by the time ordered exponential 

U(t) ^Ve^ A(t ' ]dt ' . (24) 

Alternatively, we can chose U(t) as an arbitrary one pa- 
rameter (differentiable) family of SL(2, R) matrices, and 
compute the Lagrange multipliers by derivation. The 
dynamics of the (if, v) sector is the same, with the same 
U(t) (one has only to remember that v appears second in 
the (if,v), unlikely u). This gives the complete solution 
of the classical equations of motion. 

In conclusion, the general solution of the Lagrange 
equations is 



with 



u(t) = a(t)u Q + b(t)p , 
v(t) = c(t)if + d(t)v . 



a(t)d{t) - b(t)c(t) = 1 



(25) 



(26) 



and (uo, ^chPcb ^o) satisfying the constraints, that is p$ 



and uo ■ Po — vo ■ t?o- The corresponding 
Lagrange multipliers are obtained from (p3|). 

N(t) = b{t)a{t) - a(i)b{t) , (27) 
M(t) = c(t)d(t) - d(t)c(t) , (28) 
X(t) = d(t)d(t) - b(t)c{t) . (29) 

As expected for a fully constrained system, a solution 
of the equations of motion is given by a one-parameter 
family of gauge transformations. 

Let us construct the general solution in a given gauge. 
We consider the gauge M = — 1, N = +1 and A = 0. 
The matrix A is then the unit antisymmetric matrix (and 
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time independent) and its holonomy U(t) is the rotation 
matrix by an angle t. We still have three arbitrary gauge 
fixings to impose at t — 0. We choose u 2 = if 2 , u ■ p = 
and it 2 (0) = 0. Using the constraints and the general 
solution (E5I), we obtain 



u(t) = (rcos(et) , rsin(et)) , 

v(t) = (r cos (et + (f>) , r sin (e't + <p)) , 

p(t) = (— resin(ei) , re cos (et)) , 

7f(t) = (re' sin (e't + (f>) , -re' cos (e't + <p)) , 



(30) 



with e = ±1 and e' = ±1. In this gauge, the two vectors 
u and v have the same length and rotate with the same 
angular speed, equal to one. Notice that the solution 
depends on two (continuous) parameters, r £ R + is the 
length of the vectors, and G Si is their relative angle 
at t = 0. Since the space of solutions is two-dimensional, 
there is a single degree of freedom, as anticipated. In 
addition, there are the two discrete parameter e and e'. 
These distinguish four branches of the space of solutions, 
in which each of the two vectors rotate either clockwise 
or anti-clockwise. 



where <?y is the diagonal matrix [1, 1, —1, —1]. From this 
observation, it easy to compute the Poisson algebra of 
the Oij observables 

{Oij, Om} = gikOji - guOjk + g 3 iO lk - gjkOu. (34) 

Therefore the Poisson algebra of the six gauge invari- 
ant observables Oy is isomorphic to the Lie algebra of 
50(2,2). 

Since the physical space is two-dimensional (1 degree 
of freedom), there are at most two independent continu- 
ous observables. Therefore there must be four relations 
between the six observables Oij , when the constraints are 
imposed. These relations can be easily obtained by com- 
puting the observables Oij in the gauge fl3(i| ) at t = 0. In 
fact, a relation between gauge invariant quantities which 
is true in a particular gauge is also true in general. From 
(pp|) we have 



Ol2 = 
013 



eJ, 

= J COS ( 

= J sin (t 



34 
O23 



e'J, 

ee'Jcos^. 
— ee'Jsint 



(35) 



where we have introduced 



III. OBSERVABLES 

Dirac Observables. An observable is a function on 
the constraint surface that is invariant under the gauge 
transformations generated by all first class constraints. 
Equivalently, an observable is a function on the phase 
space which has weakly vanishing Poisson brackets with 
the first class constraints. To find gauge invariant ob- 
servables, we can proceed as follows. As already no- 
ticed, ( |l4| ) indicates that the four two-dimensional yec- 
tors Xi = (xl^xf), i = 1,2,3,4 



P 



J'2 



P 











), £3 = 1 




,£4 = 1 


W 2 ) 



transform under gauge transformation in the fundamen- 
tal representation of SL(2,R). But SL(2,R) preserves 
areas in R 2 , that is, it preserves the vector product of 
any two vectors. It follows immediately that the six ob- 
servables 



Oij - — X Xj — X^ X j ) 

are all gauge invariant. Explicitly: 



(31) 



012 


= «y 


-p'u 2 , 


023 


= uV 


-A 1 , 




013 


= u'v 1 


-A 1 , 


024 


= uV 


-A 2 , 




Ol4 


1 2 

— U V 


-A 2 , 


034 


= 7rV 




(32) 



The Poisson brackets between the components of the x% 
are 

{xj,x)} = 0, {x 2 ,x 2 A = 0, {x\,x 2 } = g ir (33) 



(36) 

Clearly 

(37) 
(38) 
(39) 
(40) 

In the last equation, indices are raised with gy. Since the 
Oij are gauge invariant, these relations hold in general on 
the constraint surface. 

Thus, the two continuous quantities J € R + ,4> G <Si 
and two discrete quantity e,e' = ±1, defined in general 
by ©, namely by 



J = 


r 2 . 


e 034 = 


e' Oi2 , 


e 2 4 = 


e' O13 , 


e 023 = 


-e' Ou 


(>,,<>" = 


0, 



1 2 

up 



1 2 



p 1 u 2 \ 



\n 1 v 2 — w 1 ^ 2 ! 



1 1 2 



= arctan ■ 



- P V| 

1 2 

u v - 



(41) 



are gauge invariant observables. They can be taken as 
coordinates of the physical gauge-invariant phase space. 
Using (^J,^5|,^0|) , straightforward algebra yields the phys- 
ical Poisson brackets: 



{J,0} = ee'. 



(42) 



(e and e' commute with everything.) Notice that J = 
is a single point (whatever 0, e and e'). Therefore the 
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phase space has the topology of four cones connected at 
their vertices (J = 0). See Figure 1. 



R : = e Je* = e(0 13 + tO M ) = e'(0 24 - *O aa ), (47) 
S := eJe-** = e(0 13 - *0 M ) = e'(0 24 + 4 23 ), (48) 




FIG. 1. The topology of the phase space. 
Notice that 

Oi2 = e J = u x p , 
-Os4= -e'J = vXTr. (43) 

are the "angular momenta" of the two two-dimensional 
"particles" u and v. Since, from @, (Cq 2 ) 2 = (-0 34 ) 2 , 
the two particles have the same "total angular momen- 
tum" . In the gauge (|30| ) , w and u rotate at equal angular 
speed: each one of the four cones represents an orienta- 
tion of the two rotations, J is their angular momentum 
and 4> determines relative angle between u and v. 

The other four 0%j arrange naturally in a 2 x 2 matrix 

M ab = ( ° 14 ^ = u a V» - p a ir b , (44) 

\U23 l^ 2 4 / 

where a,b = 1,2. If we solve ( p5| ) for a(t) , b(t) 7 c(t) and 
d(t) and we insert the solution in (p6[), we obtain with 
some straightforward algebra 

u a (i) V b (t) e ac e b(i M cd = 12 34 . (45) 

(The Ojj and M cd observables are time independent.) 
Using (pq), this relation becomes 

[V^u 1 ^) + eeV(i)v 2 (i)] cos0 
+ [V(t)t; 2 (i) - eeV^u 1 ^)] sin0 = J. (46) 

This is a key equation, which entirely captures the phys- 
ical content of the model. It expresses the relation be- 
tween the Lagrangian variables (u,v) in each (J, </), e, e') 
state. The state of the system, (J, 0, e,e'), cannot be 
computed from the knowledge of the position u, v at a sin- 
gle time: two times, or a time derivative, are needed, as 
for any dinamical system. Once the state is determined, 
equation (^) provides us with the entire gauge invariant 
information: the relation between the Lagrangian vari- 
ables at any other time. 

We also define the two complex conjugate observables 



which will be convenient in the quantum theory. A com- 
plete set of observables is given by J, R, S, e, e' with the 
reality conditions 

1 = J, R=S. (49) 

Clearly 

cos (fi = —(R + S) J -1 , sin0 = —(R — S) J -1 . 
2e 2ie 

(50) 

Evolving constants. The physical phase space is the 
two-dimensional space of the gauge orbits on the con- 
straint surface. A point in the physical phase space is 
determined by (J, <fi, e, e'). This description of the system 
resolves gauge invariance, but looses reference to time 
evolution. Time evolution is, as in any fully constrained 
theory, a gauge transformation. 

In certain fully constrained physical models such as the 
free relativistic particle or the Nambu string, there is a 
global implementation of the kinematical Poincare group. 
The generator of this group that corresponds to the en- 
ergy, can be taken as the physical Hamiltonian for time 
evolution. In other words, for these systems the natural 
time evolution can be introduced in the frozen reduced 
phase space by using the energy as Hamiltonian. This 
provides a preferred variable that plays the role of time, 
namely of the independent evolution parameter. Instead, 
the kinematical group is absent in GR (unless additional 
structure, such as flat asymptotic infinity is added), or 
in the model studied in this paper. In these cases, there 
is no preferred time variable. The theory just describes 
-very democratically !- the relative evolution of the vari- 
ables, as functions of each other, without privileging any 
variable as the independent one. For a detailed discussion 
of the physical meaning of this very important feature of 
GR, see @. 

One way to express evolution in these cases, is to break 
gauge invariance. For instance, one can impose a time 
dependent gauge fixing (the analog of x° — t for a rela- 
tivistic particle), or choose a gauge at time zero and then 
evolve with arbitrarily fixed Lagrange multipliers. This 
amounts to arbitrarily choosing one of the variables as 
the time variable. 

Is there, in alternative, a gauge invariant description 
of time evolution? Are there gauge invariant observables 
that capture the dynamics of the Lagrangian variables 
u(t),v(t)7 Can we talk about a gauge-invariant dynam- 
ics, if the time dependence of u(t) and v(t) is a gauge 
transformation? The answer is yes || . 

In fact, the gauge invariant (or physical) content of the 
model is not the description of the evolution of the 4 real 
variables M 1 (t), u 2 (t), v 1 (t), v 2 (t) in the coordinate time t, 
but rather the description of their evolution as functions 
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of each other. More precisely, since there are 4 variables 
and the gauge orbits are 3-dimensional, the system de- 
scribes the motion of any one of these four variables as 
function of the other three. In other words, once the state 
of the system is known, the dynamical model allows us to 
predict the value of any one of the four Lagrangian vari- 
ables from the value of the other three. This prediction 
is uni vocal and gauge-invariant.^] 

Each solution of the classical system, namely each 
point of the phase space determines one functional re- 
lation between the four variables u 1 (t) , u 2 (t) , v 1 (t) , v 2 (t) . 
This functional relation allows us to compute one of these 
variables from the value of the other three. This func- 
tional relation is given by equation (fifil). 

The form of a gauge invariant observable describing 
evolution is therefore the following. Let us ask what 
is the value U 1 of the observable u , when u 2 and v 



have assigned values u 



x, v 



y and v 



In other words, we search an observable of the form 
U 1 = U 1 (x, y, z; J, (j>, e, e'). Solving ( f46| ) for u 1 , and re- 
placing u 2 ,v l and v 2 with x,y and z, we obtain 



U x (x,y,z; J,(j),e,e') = 



—e'x(z cos 4> — y sin cf) + eJ 
e(ycos0 + zsin <p) 



(51) 



This is an "evolving constant" in the sense of refer- 
ence [0. For any fixed state (J,<p,e,e'), the quantity 
U (x,y,z;J,(fi,e,e'), viewed as a function of x, y and 
z gives the evolution of u 1 as a function of the other 
variables. Viceversa, for any fixed x, y, z, the quantity 
U 1 (x, y, z; J, 4>, e, e'), viewed as a function of J, 4>, e and e', 
defines a gauge invariant function on the physical phase 
space. Similar expressions can be derived from (|4^) for 
u 2 , v 1 and v 2 . 

U 2 {s, y, z; J, <j), e, e') 
V 1 {s,x,z; J,(p,e,e') 

ts cos qj — t jo am ip 

2 , -y{es cos cj> - e'xsincj)) + eJ 
V (s, x, y; J, 0, e, e )= — — , (52) 



— es(y cos (j) 


+ z sin 4>) + 


eJ 


e'(zcosi 


p — y sin (p) 




— z(e'xcos<j. 


> + es sin <f>) 


+ eJ 


es cos (, 


!> — e'x sin </> 




— y(es cos cp 


— e'x sin tfi) 


+ e,J 



e x cos < 



es sin < 



where s is the value of u . These observables describe 
the evolution of the system and are gauge invariant. 

Time reparametrization invariance. The system is in- 
variant under time reparametrization. If (u(t),v(t)) is a 
solution of the equations of motion, then 



u'it) 
v'(t) 



v(f(t)) 



(53) 



*The situation is exactly the same as in GR, where the the- 
ory does not allow us to predict the value of the fields at given 
coordinate values, or the coordinate positions of particles, but 
rather allows us to predict the value of general covariant quan- 
tities, such as the value of the fields when (and where) certain 
other dynamical variables have given values M. 



is also a solution. This is immediately seen from (|25|) and 
©, because a(f(t))d(f(t)) - b(f(t))c(f(t)) = 1 follows 
from a(t)d(t) - b(t)c(t) = 1. 

Notice that there exist gauges in which u(t) evolves in 
t while v(t) remains constant. For instance we can choose 
M = A = 0. In this gauge, 



1 = ' o T) 



and therefore 



so that 



U(t) 



b(t) 
1 



u = u + b{t)p Q , 

v = v 



(54) 



(55) 



(56) 



with N — b. A different example is the following. The 
solution (^0|) can be gauge transformed to the solution 

v = (u , 0) , 

n = (u sin£ , 0) , (57) 



u = (u cos t , u sin t) , 

^ / (cosi-1) 

P = U — , u 

\ smi 



where the Lagrange multipliers are A = co !* , M — 

— ^7^^j2 , and N = 1. Similarly, there is a gauge in 

which v(t) evolves in t while u(t) remains constant. 

Notice, however, that there isn't really a "two finger 
time reparametrization invariance" in the system [^|, in 
the sense that it is not true that if (u(t),v(t)) is a solution 
of the equations of motion, then 



u'(t) 
v'(t) 



(58) 



is also a solution. In fact, in any given time u'{t) and v'(i) 
must be connected to the same point in phase space by 
a gauge transformation, but in general it is not true that 
a{h{t))d{f 2 {t)) - b(h(t))c(f 2 (t)) = 1 when a(t)d(t) - 
b(t)c(t) = 1. 



IV. QUANTUM DYNAMICS 

We work in the coordinate representation. Elements of 
the Hubert space are functions ^(u, v) of the coordinates, 
and the momentum operators are 



P : 



-ih\7 u , TT — ~iTiS7 v . 



(59) 



By inserting these operators in the constraints we obtain 
the Dirac quantum constraints 



Hi 



1 



(h 2 A v 



H 2 = --(h 2 A v +u 2 



D 



-in ( u ■ V u — v • V, 



(60) 
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where A u = V?, = 



d 2 



In the Hamiltonian 



constraint operators H\ and H2 there is a natural order- 
ing. In the "diffeomorphism" operator D, we have chosen 
the ordering that leads to the closure of the constraint 
algebra and thus the absence of anomalies. We have in 
fact 



H\ , Hi 
H X ,D 
H 2 ,D 



= ihD, 
= -2ihH 1 , 
= 2ihH 2 . 



(61) 



The physical states, in the sense of Dirac, are in the 
kernel of all the quantum constraints. Namely, they are 
defined by 



v ) "i>(u,v) 



0, 



(n 2 A u 

(h 2 A v +u 2 ) V(u,v)=0, 
-ih (u-'Vu-v- V v ) v) = . 



(62) 



We now solve this system of coupled partial differential 
equations. 

We transform to polar coordinates 

5= (ucosa, usin a), v = (v cos/3, vsin/3) (63) 

and we multiplying the first equation of the system by 
u 2 /h 2 and the second by v 2 /fi 2 . becomes 



d_ d_ 

du du 

d_ d_ 

dv dv 



d 2 
da 2 

d/3 2 ' 
d_ 
du 



dv 



*(w, v,a,/3) 
ty(u, v, a, (3) 



(64) 



We search a solution by separation of variables, by writ- 
ing 

*(tt, v, a, 0) = A(a) B{f3) tp{u, v). (65) 

The first two equations in (^34j) give immediately 

A(a) = e™° Q , B(fl) = e™ 3/3 , (66) 

where m a and mp must be integer for \I/ to be continuous. 
The third equation in (t34j) implies that 



ip(u, v) — ip(uv) 



(67) 



(a function of the product uv). Plugging this last result 
back into the first two equations in (p4j), we find that the 
first and last terms of one equation are equal to the first 
and last terms of the second. Therefore the two middle 
terms must be equal as well. Therefore the two equations 



m a = em, m/3 



-em, e, e 



±1. 



(68) 



with m any nonnegative integer. The minus is inserted 
for later convenience. Using this, the first two equations 
of the system become equal to each other and reduce to 



d 2 f{x) 
dx 2 



1 df(x) 
x dx 



1 - 



/(s)=0, 



(69) 



where we have written x = uv/h and f(x) = ip(hx). 
This is the Bessel equation. Thus, we have solved the 
system entirely. We conclude that the physical Hilbert 
space is spanned by the basis states \m, e, e'), where m is 
a nonnegative integer and e, e' = ±1. In the coordinate 
representation these states are given by 



(u,v,a,(3\m,e,e') = 9 m ^ e >(u,v,a,j3) 
- e J m y h 



(70) 



where J m is the Bessel function of order m. Notice that 
for each m > there are four states (e = ±l,e' = ±1), 
but for m = there is only one state, since |m, +, +) = 
\m,+,-) = \m,-,+) = \m,-,-) »( 

Quantum observables and scalar product. Consider the 
observables Oy defined in (32). They are gauge invari- 
ant, and thus have vanishing Poisson brackets with the 
constraints. We chose the natural ordering for the corre- 
sponding quantum operators Oy 



6 12 


~ 1 ~2 

= up 


A 1 ~2 

— jrti , 


6 23 


= uH l 




Ol3 








= u 2 v 2 


-P 2 K 2 


Ol4 


= uH 2 




034 




-v l % 2 



(71) 

It is easy to see that the commutators of these opera- 
tors with the quantum constraints (60) vanish. Therefore 
these operators are well defined on the space of the so- 
lutions of the quantum constraints, namely on the states 
(f70|). We compute their action on these states. Going to 
polar coordinates we see immediately that 



d 

O12 * m , e ,6' = -ih—% 

da 



O34 * 



34 ^m,e,e' 



(72) 



Thus in the physical state space we have e'0\2 = eO^: 
the relation between the two is precisely the same as in 
the classical theory, eq. (37). We can thus identify the e 
and e' appeared in the quantum theory with the e and 
e' appeared in solving the classical theory, and we con- 
clude, from equation (35), that the quantum operator 
corresponding to the gauge invariant observable J is 



imply m c 



m r . 



We put 



'We missed this point in the first version of this paper. We 
thank Jorma Louko for pointing out the mistake. 
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J \m, e, e ) = hm |m, e, e ) 



(73) 



Thus in the quantum theory J is discrete, quantized in 
multiples of h 

J = mh. (74) 
Using the Bessel equation and the properties 

x ax 



X ax 



(75) 



of the Bessel functions, a straightforward but long calcu- 
lation yields 



(O13 + iOu) ^m,e,e> = ehm vt/ m+ee , e e , 

(0 2 4 - i0 2 3) *m,e,e' = e'^77l ^m+ee'^^ 



(76) 



Thus, the quantum operator corresponding to the ob- 
servable R defined in ( (47| ) is 

R |m, e, e ) = fi.rn \m + ee', e, e'). (77) 

In the same manner, from ( (48| ) we obtain 

5 |m, e, e') = Tim \m — ee', e, e'). (78) 

To complete the construction of the Hilbert space of the 
physical quantum states, we have to determine the scalar 
product on the space spanned by the states \m, e, e'). 
This is determined by the requirement that real classi- 
cal observables be self adjoint. The observables J, e and 
e' are real, and thus we require J, e and e' to be self ad- 
joint. It follows that the states |m, e, e') which are their 
eigenstates must form an orthogonal basis. This fixes the 
scalar product up to the norm of the basis states. Define 



(m,e, e'\m, e,e') = c n 



(79) 



Next, S is the complex conjugate of R. Thus we require 
that R f = S. It follows 

(m, e, e'| R^ \n, e, e') = (m, e, e'| S \n, e, e'). (80) 

From which, we have easily 

Cm.e.e' = Cm. (81) 

Here c is a positive overall normalization constant that 
has no effect on the physics, and we chose equal to 1. This 
fixes the normalization of the orthogonal basis states, and 
therefore determines the scalar product completely. No- 
tice that the state |0, e, e') has zero norm. (This was first 
realized by Jorma Louko). We can therefore discard it, 
because its presence has no physical consequences. More 
precisely, we identify the m = state with the state zero. 

The peculiar behavior of the m — sector of the quan- 
tum theory reflects the pathological properties of the cor- 
responding classical state. The quantum state m = has 



vanishing angular momentum J; the classical state with 
vanishing angular momentum is the (common) vertex of 
the four cones that form the reduced phase space (see 
Figure 1). This is a point at which the reduced phase 
space fails to be a manifold. Physically, this corresponds 
to the fact that small perturbations of the J = solution 
form disjoint spaces. 

Thus, the quantum theory is completely defined by the 
states 



m — l,oc;e,e'— ± 



|m,e,e'), 



the scalar product 

(m, e, e'|m, e, e') = m (5 mi m^e,e^e 
and the operators 



(82) 



(83) 



J 


m 


e, 


e ) 


= hm \m, 






R 


m 


e, 


e ) 


= hm \m 


+ 


ee', e, e') 


S 


\m 


e, 


e ) 


= hm \m 




ee',e,e') 


€ 


m 


e, 


e ) 


= e \m, e 


e 


>, 


e' 


\m 


e, 


e ) 


= e \m, e 


e 


)■ 



(84) 



where it is understood that |0,e, e') = 0. (That is, for 
instance, R |1, +, — ) = 0). 

Quantum evolving constants. In order to quantize the 
evolving constant of motion ([51 52), we must construct 
the operators corresponding to the classical observables 
cos<^ and sinc/i. We denote these operators cos0 and 
sin </>, with a slight abuse in notation. (The operator 4> 
is ill defined because <f> is an angle -see for instance ||- 
and we must deal with periodic functions of (j) in order 
to have continuity all around the circle.) Choosing the 
natural ordering given in (|50"|), we have immediately 

cos0 |m, e, e') = 7p(l m + e , e') + \m— ee', e, e')), 



sm0 m, e, e 



1 

2d' 



m + ee , e, e 



')-|m-ee',e,e')) (85) 



(where, again, it is understood that |0,e,e') = 0.) 

A convenient representation of the theory can be ob- 
tained by representing a generic state 



\i>) = ^2 \m,e,e') 

m,e,e' 

by the four functions on 5i 

oo 

The scalar product turns out to be 

f d 



(86) 



(87) 
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Notice that since the sum in (J87j) is restricted to m > 0, 
the Hilbert space is formed by "right moving" func- 
tions ?/>+,+ (</>) and ip- t -((f>), and "left moving" functions 
ip+ t —(<f>) and only. On these functions, the scalar 

product (§8|) is positive definite. In particular, the zero 
modes ip e ,e' (0) = constant do not belong to the Hilbert 
space. We denote the projector that projects out the zero 
modes as P. The observables are then 



J Tpe,e'(<f>) = 



-ihee' -Tjtpe 



COS0 1p e , e /((f)) = P COS0 1p eye r((f)), 
Sin0 4> e ,e'{4>) = P Sin*/ 1 4>e,e'{4>)- 



(89) 



In this representation it is easy to write the quantum op- 
erator corresponding to the evolving constant of motion, 
which quantizes the observable d5l|). This is given by 



P- 



-ee 



y cos <p + z sin < 



-P 



x(z cos 4> — y sin 0) + ih— 



where we have arbitrarily picked an ordering.^ The ex- 
pectation value of this operator on a state ^^^{(j)) - 
taken with the scalar product (p8|)- gives the physical 
mean value of the variable u 1 at the moment in which 
the three variables u 2 ,?; 1 and v 2 have value x,y and z 
(see ||). Similar operators can be defined for the three 
other evolving constants 
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